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We present a completely gauge-invariant operator definition of transverse-momcntum- 
dependent parton densities (TMD), supplied with longitudinal lightlike gauge links as 
well as transverse gauge links at lightcone infinity. Within this framework, we consider 
the consistent treatment of specific divergences, emerging in the "unsubtracted" TMD 
beyond the tree approximation, and construct the soft factors to cancel unphysical sin- 
gularities. We confront this approach with factorization schemes, which make use of 
covariant gauges with off-the-lightcone gauge links, and discuss their mutual connection. 

1. Introduction 

Different operator definitions of the transverse-momentu m-dependen t parton dis- 
tribution functions (PDF) — TMD for short in what -are actively 
discussed in the literature, see, e.g., Refs. [3 El El [TU1 HH HH H31 HS1 [TH1 H3 [TBI 
and references therein. Among the principal issues to be addressed in any consis- 
tent operator definition of the TMD, there are their gauge invariance, the cance- 
lation of undesirable divergences, renormalization-group properties, and ultraviolet 
(UV) and rapidity evolution equations. Moreover, properly defined TMDs have to 
be incorporated into the QCD factorization formula for the structure functions of 
semi-inclusive processes that can be schematically written as 

F S emi-incl. ~ H (g> T D ® T F ® S , (1) 

where H is the hard (perturbatively calculable) part, J-d,f are the TMD distribu- 
tion and/or fragmentation functions, and S is the soft factor — a specific ingredient 
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in the semi-inclusive factorization approach. Nonperturbative distribution functions 
of partons (in what follows we consider only quark distributions) , depending on the 
longitudinal components x, as well as on the transverse components k±, of their 
momenta accumulate information about the intrinsic motion of the hadron's con- 
stituents. Initially, the TMDs have been considered as a direct generalization of the 
collinear PDFs: 



F[ co n.](x, n) -> F[ tmd ](x,^±,/j,,ri) , / d k±F [tmd] (x,'k±,ti,'n) -)■ F [colL] (x, fi) . (2) 

However, such a straightforward relationship between TMDs and integrated PDFs 
can only be proved in the tree approximation because rapidity divergences jeopar- 
dize this procedure, or even entail its breakdown. Note that in the collinear case 
the only scale is the UV renormalization parameter ji, whil e in the TMD case, an 
additional dependence from the rapidity cutoff £ arisei^EEl 

In this work, we present and discuss an operator definition of the quark TMD 
that embodies gauge invariance in terms of lightlike longitudinal and transverse 
gauge links with the appropriate behavior at lightcone infinity ^ 



2. Gauge links in the lightcone gauge 

We start with the "unsubtracted" definition (i.e., without the explicit isolation of 
the soft factor) of the "quark in a quark" TMD (where the subscript on [A] n denotes 
the axial lightcone gauge for lightlike longitudinal gauge links along the vector n, 
employing the notations in Ref. 20): 

1 fdrd 2 ti± 



it„,+ 



27r(27r) : 



--e- ik < (pi&(r,€-L)[r,ej-;«>-,£i]t 



x[oo ,£x;°o ,oo ± ]|7 + [oo ,oo ± ;oo ,0j_]j[oo ,0±:0 ,0x] n ipi(0 ,0±)\pJ(3) 

with £ + = 0. The path-ordered longitudinal (lightlike, n 2 = 0) [...]„ and transverse 
gauge links 



[oo ,^_L]n = ^exp -ig 



oo_l;oo = Vexp 



dr n~ A£r(t + TTT) 

POO 

-ig / dr I- A a t a (£±+lT) 
Jo 



(4) 



ensure the formal gauge invariance of this TMD. The TMD can be normalized as 
follows: 

1 f rfr^X -ifc+£-+ife,.£, 

27r(27r) 2 



r [A„](0) 

unsubtr. 



x W(r,£±)7 + ^((r,0 ± )|p) = 6(1 - x)S^(k x ) 



(5) 



In the tree approximation (which we distinguish from ^°„ b ° t i. by keeping "clas- 
sical" gauge links in the former) straightforward integration over the transverse 
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momentum fej_ immediately yields the collinear gauge-invariant PDF 

= \J ^t e ~ lk+r <^i(r,Ox)[r,0-]„ 7 +^(0- Ox)b) . (6) 

Let us emphasize that the ofF-the-light-cone operator definitions (A v -, G V -TMD 
given in Ref. I20() do not satisfy this relation — even in the tree approximation. In- 
stead, they can be shown to obey a factorized formula in terms of collinear PDFs 
at small b± in the impact-parameter representational. The use of one of these 
approaches is a matter of convenience, provided that the corresponding operator 
definition of the TMD is consistent with the factorization scheme, cf. Eq. ([I}. How- 
ever, one has to be careful when comparing different classes of definitions, like A v -, 
C V -TMD vs. A n -, C n -TMD, since these are, in principle, different objects even at 
the tree level. In contrast, the comparison of the definitions A n vs. C n is justified, 
since they represent the same quantities in different gauges. 

Beyond the tree level, the quark fields in the operator definition of the TMD ^ 
have to be considered as Heisenberg field operators, i.e., 



dr\ ipAip 



A ^(vhMv)^(v) • (7) 



Therefore, the O(g) contributions from the gauge links are contracted with the 
quark-gluon interaction terms J drj ipAip , originating from the Heisenberg fields 

(O, and give rise to the set of the 0(g 2 ) one-gluon exchange graphs shown in Fig.Q] 
(more details are given in Ref. |2T|) . 

In the light-cone gauge, some of these diagrams disappear. Let us show that 
the contributions of the longitudinal (lightlike) gauge links (which are shown in 
Fig. Q] by horizontal double lines, whereas the transverse gauge links are denoted 
by vertical double-lines) cancel in the lightcone gauge A + = when combined with 
^-dependent pole prescriptions for the gluon propagator. Formally, the "classical" 
exponential exp[— ig Jdz~ A + ] equals unity in the lightcone gauge, irrespective of 
the applied pole prescription, but this is worth being proven for the prescriptions 
we use in our analysis. 

Next, we evaluate the gauge field, the source of which is a charged pointlike 
particle (a struck quark) moving with the quasi-constant four-velocity along the 
straight line v^t. The corresponding "classical" current then is 

Mv) = 9 fdy' s {4) (y -y'), y'u = ^r. (8) 



The velocity changes its direction only at the origin, where the collision with the 
hard photon takes place and the quark deviates from its initial "trajectory". The 
gauge field defined by such a current then reads 

^ c (0 = Jd*yv&(Z-y)My) . (9) 



4 Cherednikov & Stefanis 




+ (h.C. 



+(h. 




Fig. 1. Complete set of the one- loop diagrams corresponding to the gauge-invariant quark TMD 
PDF without soft- term contributions. 



where is the free gluon propagator in the lightcone gauge. This gauge field is 
exactly what forms the longitudinal gauge hniia We assume that the velocity of 
the struck quark is parallel to the "plus"- and the "minus"- lightcone vectors 
before and after the hard collision, respectively: 



J/j,(y) = 9 



Then, one has 



n+ I dr 5^ (y - n+r) + n~ / dr 6^ (y - n+r) 
dq~ e-*~* 



K s (y ) 



2tt q- 



iO 



n 8{y^ 



dq + e~ iq+v 
~2n q+ - iO 



.(10) 



~9 nt 



d A q 



e-«* V&(q) dy + dy-d 2 y ± e^yS(y-)S^(y ± ) , (11) 



2(2tt) 4 " " LCV 

where the free gluon propagator in the lightcone gauge A + = (A ■ n~) = reads 



d 4 q 



,q - z V^(q) = - 



d A q e 
(2^)4 If 



it) 



q^{n-y + q v (n-y 



(12) 

and the pole prescription [q + ] v has yet to be defined. We neglect the quark and 
gluon masses, since we are mainly interested in the UV and rapidity singularities. 
After the integration over the variable y, we get 
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Now we are able to calculate the plus-component of the gauge field (|TU| A + 
(A-n-y. 



1 f d 1 „-ia-£ % ) 



(14) 



(2tt) 4 g 2 + iO 

The first integral is trivial and yields (using dimensional regularization w = 4 — 2e 
for the integration over transverse degrees of freedom) 

d^q __ iq< 8{q-) 1 f d u - 2 q e iq ^^ 



(2tt)" 6 ^ c/ 2 + zO 2tt^ ) 7 (2tt)-- 2 -q\ + iO ' (15) 
To perform the ^-dependent integral, we observe that 

d u q e-*K q+ .8 f d"q e^< 5{q-) 



(2ir) u 9 2 + iO [g+]„ vy ; J (2^) w <7 2 + «0 [q+] v 

Making use of 5(q~) in Eq. (|16|) . we separate out the transverse part 

d w g e- l « « 1 r , _. 1 fdq+ f d"- 2 q e"^^ 



(2tt) w 9 2 +i0 [q+] v KH ' 2tt y 2tt [q+] n J (2tt)"- 2 -q 2 ^ + iO 
Consider now the longitudinal integral which we define as 

lim / . (18) 



27r [#+],, tj— >-q 7 27r q+ ±ir) 

where ± corresponds, respectively, to the retarded and advanced prescription, while 
the principal-value prescription can be obtained by symmetrization. Integral (|18[) 
can be evaluated using the residue theorem to read 

[^ e -=^=Ti6(±n^-\ e {z ) = {^ Z> l . (19) 
y 2tt q+ ± if] ^ V ^ ; ' W \ , z < V ; 

Taking the derivative with respect to £ _ and employing Eq. (|16j) . we find in the 
limit 77 — >■ the expression 

._d_ f d"q e-W S(q-) _ 1 f d"- 2 q e^ ^ 

y (27r)«- 9 2 + z0 [q+] v 2tt ° U V (2tt)"- 2 q\ + zO ' 1 j 

Reversing the sign and adding this result to Eq. (IT51) . we get, according to Eq. (|14l) . 

A+ = (A-n - ) = 0. (21) 
The "minus" -component of the gauge field 

A- - ■ <) ~ < ((„+)> - - C - 2 f } (22) 

vanishes as well after carrying out the integration over q~ in Eq. ()13j) by virtue of 
6(q~), i.e., A" = 0. 

This suffices to show that the use of the ^-dependent pole prescription in the 
gluon propagator is consistent with the lightcone gauge A + = 0, while for a further 
discussion we refer to Refs.[Hl[22l A formal proof of a similar statement concerning 
the Mandelstam-Leibbrandt pole prescription will be presented separately. 
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3. Lightcone TMD: One-loop effects 



In the previous section, we have shown that in the lightcone gauge with 77-dependcnt 
pole prescriptions for the gluon propagator, the longitudinal (lightlike) gauge 
links are equal to unity: exp[— ig Jdz~A + ] = 1. Therefore, the diagrams Fig. [1] 
(b, d, e, /, h,j, k) (and their Hermitian conjugate parts — if any) give zero contribu- 
tions in the lightcone gauge, despite the opposite claims by Collins in Ref. [TU 
Moreover, we give below a formal justification of our statements by demonstrating 
the gauge invariance of our framework. 

In the lightcone gauge, the 77-dependent contribution of diagram 1(a) is 



E^fAn 



C F F(e) 



4tt 



S(l-x)S^(k ± ) 



(23) 



In a covariant gauge (definition C n ), the counterpart of this contribution stems from 
the vertex diagram 1(6) and reads 



s (c) [c„] = --c F r( £ ) 



4-7T 



2 1 e 



~P 



6(l-x)6W(ki_) j^dx {1 ~_ x)1+e 



(24) 



By a trivial change of variables and by using the ry-regularization, the singular 
integral above can be written as 



£ (c) [c n ] = -^c F r( e ) 



4tt 



6(1 -x)5 (2) (k ± ) 



x £ (x±irj) 



x) 1 



-£ (a) [A n ] 



(25) 

This is an interesting result and supports the validity of the operator definition of 
the TMD with the lightlike longitudinal gauge links, justifying its gauge invariance 
at the one-loop level. Analogous results have be en obtain ed recently within the 
framework of soft collinear effective theory (SCET p^4|25 | Similar problems have 
been addressed and resolved long ago in Ref. HSj 

Therefore, whatever gauge is adopted, the one-loop corrections to the "unsub- 
tracted" definition ([3]) give rise to pathological overlapping divergences that com- 
prise UV and rapidity poles simultaneously ~ ^ In rj. This jeopardizes the renor- 
malizability of TMDs and calls for a certain generalized renormalization procedure 
to augment the insufficient dimensional regularization. In our workjA^, we worked 
out such a procedure that enabled us to obtain a well-defined and fully gauge- 
invariant TMD PDF, free of undesirable divergences. This TMD PDF has calculable 
renormalization-group properties and obeys an evolution equation with respect to 
rapidity in the impact parameter space. Moreover, the one-loop analysis of the UV 
anomalous dimension of the "unsubtracted" TMD PDF in © with the ^-dependent 
pole prescription in the lightcone gauge shows that the overlapping singularities pro- 
duces a correction to the anomalous dimension that can be identified at this order 
with the well-known cusp anomalous dimensior9^. The upshot of this discussion 
is that in order to renormalize the "naive" TMD PDF, given by (|3]), one has to 
introduce an additional renormalization factor, which depends on fj = i]/p + and 
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can be written as the vacuum average of the gauge links evaluated along a special 
contour with an obstruction (cusp) — see Ref . [TTJ 032 for more details: 



z- 1 = (o 


Vexp 


ig \ 









dC t a A a JQ 







(26) 



However, the overlapping singularities are not the only ones that have to be 
removed from the proper definition of the TMD. The one-loop corrections to the 
soft factor itself give rise to another type of unphysical divergences. For instance, 
in the lightcone gauge, the self-energy contribution of the soft factor still contains 
an uncanceled singularity of the term 

d u q 1 



J soft 



[An] 



ig 2 /i 2£ C F 2p 4 



iO)[q+] 



(27) 



(2tt) w q 2 {q- 'p-t 

which appears to be rapidity-independenlP^l. This calls for an additional subtraction 
of this self-energy part that is presented graphically in Fig.[2J This subtraction does 
not affect the rapidity evolution equations and does not break the factorization 
structure fl}. In fact, it has an intuitively clear physical interpretation: it serves 
to remove unobservable contributions due to the self-energy of the infinite lightlikc 
gauge links that are mere artifacts of the unobservable background. 





(0-,oo+ ± ) 



(0-,-oo + , Ox) 



(0- -oo+ Ox) 



Fig. 2. Subtraction of the infinite self-energy contribution of the lightlike gauge links in the soft 
factor. 



The definition of the A n TMD PDF, which embodies the above requirements, 
readJlSl 

'A (*.kx,M,»W- 2 J 27r(27r)2 e 
'/i|V'i(C _ ,^x)K _ >€x;oo _ )Cx]n[° 0_ >Cj-;oo _ ,oox]] 7 + • 

x[oo-,oo_ L ;oo-,0x] i [oo-,0x;0- ! 0±] n ^(0-,0±)|/i>-i?- 1 , (28) 
with the soft factor 

R~ 1 M = 



(0| 7? exp 




■ V 1 exp 




10) 


(0| Pexp 




• V~ x exp 


.-fc-kw^C" A^ + 0_ 


10) 



where A = t a A a and the contours for the soft factor are displayed in Fig. [2] 



8 Cherednikov & Stefanis 



4. Light-cone TMD: Evolution equations 



After subtracting the soft factor, the operator definition ([28)1 is multiplicatively 
rcnormalizable and obeys the following one-loop evolution equation with respect to 
the UV scale [i: 

^ ^ (X ' k±; ^ V) = (7LC + • F[AnI (X ' k±; M ' V) ' (29) 

where 7lc — i^t^f is the anomalous dimension of the bilocal quark operator in 
the lightcone gauge. It is interesting to note that, in the lightcone gauge with the 
Mandelstam-Leibbrandt prescription, the "unsubtracted" TMD ([3]) has an anoma- 
lous dimension that is free of any undesirable contributions and reads 



{An 

unsubtr 



R~ 



d_ [a™ l ] 



(7lc + 0(a 2 s )) 



[A? 



unsubtr 



(7LC+0(^))^ubt 



(30) 



In the one-loop approximation, these results are in agreement with those given in 
Ref.El 

In order to study the rapidity evolution, we have to concentrate on the specific 
TMD singularities which only depend on the additional rapidity parameter rj, but 
else do not violate the renormalizability of the TMDs. These si ngul ar terms have to 
be resummed by means of an equation of the Collins-Soper type™. Our framework 
indeed allows the development of such a resummation procedure of the rapidity 
divergences as we now show. 

To this end, recall that the ^-dependent part ([23)1 of the quark self-energy dia- 
gram corresponding to Fig. 1(a) has the form 



sW[AJ=-^C F r(e) 



4tt 



-p 2 



S(l-x)S^(k±) 



-1 - In 



^- £ (V^-iln 2 ^ 
p+ V 3 2 p+ 



+ 0(e 2 ) 



(31) 



where the quadratic term ~ In 77 is responsible for the "rapidity" evolution. The 
linear terms stem from the remaining virtual and real graphs of the "unsubtracted" 
TMD and the soft factor. 

In order to confront our framework with the Collins-Soper rapidity evolution 
approach, we make use of the impact representation of the TMD 



J 7 (x,b ± ; t i,r 1 ) = J d 2 kjL e lb± k± F (x, k x ; M, r?) 



(32) 



Then, the Collins-Soper rapidity evolution equation (which holds for the off-the- 
light-cone A V -TMD and C V -TMD) becomes 

d 

(q£ ^[a,c v ] [x,^±;h,0 = [K v (fj,,bx) + G v (fx,Q] J[a,c v ] (x,kx;n,Q > ( 33 ) 
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where the functions K and G have the following renormalization-group properties: 

^^t^ v = _/i (^ Gv = 7cusp ' ^ 

In our approach, the analogue of the Collins-Soper rapidity cutoff is given by the 
new variable 9 = (p ■ n)jr\. Therefore, the corresponding evolution equation takes 
the form 

9^ T [An] (x, k ± ; /i, 9) = [K n {fx, b ± ) + G n (p, 0)} T [An] (x, k ± ;/x, 9) , (35) 

noting that the limit 77 — > corresponds to the limit £ — > 00 in the Collins-Soper 
approach. Here the sum K n + G n for the A n -TMD can be evaluated perturbatively 
in the small-b_i_ region. The relation fl34[) for the A n -TMD has been verified at the 
one-loop level in our previous workJC2. Explicit results for the rapidity evolution 
kernel [K n ([i, b±) + G n (fj,, 9)) will be reported in the future. 

The bottom line is: the generalized definition of the TMD PDF (28) allows one 
to derive its renormalization-group evolution with respect to the UV scale \i and 
to resum the large rapidity logarithms by means of the evolution equation (|35[) . 
formally akin to that in the Collins-Soper procedure, designed for off-the-light-cone 
quantities. 



5. Conclusions and discussion 

To conclude, we have shown that the completely gauge-invariant operator definition 
of the transverse-momentum-dependent parton densities, which makes use of longi- 
tudinal lightlike gauge links and also of appropriate transverse gauge links at light- 
cone infinity, can be consistently formulated beyond the tree-level approximation — 
at least in the one-loop order. The advantages of the presented approach are: 

• In the lightcone axial gauge, the longitudinal gauge links cancel, while the 
transverse gauge links can be eliminated by proper boundary conditions 
for the gauge fields at lightcone infinity. This provides a natural framework 
for QCD factorization schemes, where the use of covariant gauges and not 
purely lightlike gauge links is not convenient. 

• A direct relationship between the (unintegrated) TMDs and the (in- 
tegrated) collinear PDFs can be established by straightforward k±- 
integration. One ultimately gets the well-known collinear and gauge- 
invariant PDF that fulfills the DGLAP evolution equation (in the one-loop 
order) . 

• The transverse gauge links at lightcone infinity (supplemented by appropri- 
ate boundary conditions) accumulate information about the initial- and/or 
final-state interactions in the axial gauges, in which the longitudinal gauge 
links disappear. The transverse gauge links naturally shrink to harmless 
constants after performing the fc^-integration without causing a breakdown 
of the standard structure of the longitudinal gauge links in the collinear 
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PDF, while in the case of the off-the-light-cone definitions this is, at least, 
not obvious. 

• Moreover, the transverse gauge links are responsible for T-odd effects in 
the lightcone gauge that makes it possible to apply this framework to the 
study of phenomenologically important quantities, such as the Sivers, Boer- 
Mulders, Collins functions, etc. 

• A complete proof of the QCD factorization within the presented scheme — 
in particular, the evaluation of the hard part in the lightcone gauge and the 
explicit demonstration of the independence of the full structure functions 
of unphysical scales is still lacking and will be treated in the future. 
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